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Abstract. In a one dimensional shallow optical lattice, in the presence of both 
cubic and quintic nonlinearity, a supcrfluid density wave, is identified in Bose-Einstein 
condensate. Interestingly, it ceases to exist, when only one of these interaction 
is operative. We predict the loss of superfluidity through a classical dynamical 
phase transition, where modulational instability leads to the loss of phase coherence. 
In certain parameter domain, the competition between lattice potential and the 
interactions, is shown to give rise to a stripe phase, where atoms are confined in 
finite domains. In pure two-body case, apart from the known supcrfluid and insulating 
phases, a density wave insulating phase is found to exist, possessing two frequency 
modulations commensurate with the lattice potential. 
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1. Introduction 

The mean field Gross-Pitaevskii equation (GP) [H E] is known to capture the ground 
state properties of the Bose-Einstein condensate (BEG) remarkably well, both in 
three and lower dimensional configurations [H H]. In the presence of an optical 
lattice O El El [7] , this picture changes drastically, when quantum fluctuations become 
significant [9]. For deep lattice wells, fluctuation driven metal- insulator transition has 
been experimentally confirmed flUi ITTl IT2[ IT3] . However, in shallow wells, the quantum 
fluctuations are expected to be small, making GP equation useful for identifying ground 
state structures. 

For discrete nonlinear Schrodinger equation, Smerzi et al. in 2002 [H], predicted a 
dynamical superfiuid insulator transition (DSIT), which is different from the fluctuation 
driven quantum phase transition. Modulational instability, occurring in the system, 
leads to the loss of phase coherence. Hence, DSIT takes place and the superfiuid phase 
transits to the insulating state. Followed by this prediction, Gataliotti et al. in 2003 |15j . 
reported the experimental observation of the destruction of interference pattern and loss 
of superfluidity through DSIT. In the same year, Adhikari demonstrated numerically the 
loss of phase coherence and superfluidity of BEG [16j with two-body interaction. 

In this paper, we demonstrate the existence of a superfiuid density wave state 
of BEG in an optical lattice, in the presence of both cubic and quintic nonlinear 
interactions. This phase, with twice the periodicity of the lattice, ceases to exist, 
when only cubic and quintic nonlinearities are present. We have considered two- 
body interaction to be repulsive and quintic interaction to be attractive for physical 
realization. The existence of stripe phase is shown, where superfiuid matter is found 
only in finite domains. We predict that DSIT occurs in this system also, where all 
the atoms transit from the superfiuid phase to an insulating phase, with analogous 
periodicity. The filling fraction in this case, varies inversely with the lattice site n and 
reaches an asymptotic constant value. It differs for odd and even lattice positions with 
respect to a given site, taking a constant value for even n. Interestingly, for odd n, it 
depends on the depth of the lattice potential, as well as the interaction strengths. In 
the presence of quintic nonlinearity alone, a sinusoidal excitation is found, which exists 
only in the superfiuid phase, of similar periodicity as that of the optical lattice. In all 
these cases, the dispersion plays a sub-dominant role. In presence of pure two-body 
interaction, the superfiuid density shows a periodicity commensurate with the lattice. 
However, a density wave insulating phase is found to exist, possessing two frequency 
modulations commensurate with the lattice potential. The dynamical phase transition 
connects the superfiuid phase with an insulating phase [I7j, which in turn is connected 
with the density wave insulating phase. The solutions are found to be marginally stable, 
as per the Vakhitov-Kolokolov (VK) criterion |l8j. 
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2. BEC with cubic and quintic nonlinearity 

The quintic nonlinearity in one dimension, corresponding to an effective three-body 
interaction, has appeared in a number of works. Muryshev et al. [19], have investigated 
the dynamical stability and the dissipative dynamics of solitons in BEC, where the 
quintic nonlinearity arises in one dimension due to the interaction between axial and 
radial degrees of freedom. In [20], the cubic GP equation is used in one dimension, 
when the coupling constant is small, whereas, quintic interaction appears in the Tonks- 
Girardeau regime. The localization of ground state of BEC in optical lattice and the 
stability of BEC with quintic nonlinearity, have been analyzed by Abdullaev et al. 
[2T| |22] . Kolomeisky et al., took a different approach and derived the GP functional 
with cubic and quintic nonlinearity for lower dimensions, using renormalization group 
method [231E1]. 

The quasi-lD GP equation in an optical lattice, in presence of both cubic and 
quintic nonlinearity, takes the form [20112211211 125] . 

#t = -\i^zz + {qM? + 92\i^\^ + v{z) - fi)^, (1) 

where, V{z) = Vocos^(z) is the optical lattice potential, gi and g2 are the appropriately 
normalized strength of the cubic and quintic nonlinearity, respectively [26]. Since, both 
cubic and quintic terms are functions of scattering length [25], we describe the effect of 
the interactions in terms of the ratio between cubic and quintic nonlinearity k = gi/g2- 
The lattice and the chemical potentials are normalized in terms of the recoil energy 
Er- The spatial co-ordinate and the wavefunction are scaled in units of wavelength of 
incident laser light and y/k respectively; here k is the wave vector. 

The ansatz solution of Eq. ([T]), ip{z,t) = \J ^{z) exp {ix{z) — iut) leads to a 
superfiuid phase (SF): 



^SF = X ( ± ^-^ COS. (2) 




where, x{z) = c tan-^ [J -'^^l'-^^^ tan(f )] , = |(1 - M _ g^) and c 



-gi/2g2+y/-Vo/g2 ^2^J' 8 



^^ydi + ^92^0- The periodicity of the density is twice that of the lattice potential, 
which indicates matter redistribution. Atoms from neighbouring sites have been 
depleted leading to a density wave behaviour. As one can see from the expression, 
the phase vanishes where the superfiuid density wave attains its maximum value. Since, 
the superfiuid is characterized by a definite phase, the number of atoms on each lattice 
site is unknown, which allows the atom to move freely and it can easily tunnel from one 
lattice site to another. 



Denoting ^ y ~v| solution exists only in the following regimes, when the 

quintic nonlinear interaction and the lattice potential have opposite signature: 

(i) If cos z > p, then 
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a) the solution exists for all 2 for p < —1, which indicates that the cubic and 

quintic nonlinearities are of opposite signature. 

b) the solution does not exist for j9 > 1. 

c) in the domain where, —1 < p < 1, the solution exists for — 7r/2 < z < Zj < 7r/2, 

where zj = cos~^ p. This domain corresponds to a stripe phase, where 
superfluid matter is found only in finite domains, 
(ii) If cos z < —p, then 

a) the solution does not exist for p > 1. 

b) the solution exists for all z for p < —1. Therefore, the cubic and quintic 

nonlinearities are of opposite signature. 

c) when —1 < p < 1, the solution exists for —tt/2 < z < —Zj < tt/2, where 

Zj = cos^^(p). This also corresponds to a stripe phase. 

The average number of atoms per lattice site, the filling fraction, is found to be 




10 20 30 40 50 60** 

Figure 1. The behaviour of the filling fraction when both cubic and quintic nonlinear 
interactions are operative, n is the lattice location with respect to a given site. The 
solid curve shows the filling fraction for the superfluid phase and the dotted curve is 
for the insulating phase with Vq = 0.002 (unit of recoil energy) and k = 5. 

Starting from a given lattice point, usf depends on the lattice location n, as depicted 
in Fig.([T]). If n is odd, the filling fraction is a function of the lattice potential and varies 
inversely with the lattice site number n. If n is even, the filling fraction takes a constant 
value depending only on the interactions. The different behaviour of the filling fraction 
at the neighbouring sites can be traced to the presence of cos{z) term in the density. 
The odd-even effect arises due to the interplay between the cubic and quintic nonlinear 
interactions, where a density wave is present. As will be seen below, for the pure cubic 
or quintic nonlinearity, this type of behaviour is absent. 
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The energy of the superfluid component is given by, 
EsF = ^\9l-3g2{2 + Vo~2fi) 




^9.V, - I , (4) 

91 + ^92Vo 



which possesses a branch cut at, Vq = where the supercurrent (J = c) vanishes. 

At this value, the phase coherence is lost and the superfluid phase transits into an 
insulator (I) phase, which is driven by a modulational instability [H]. The insulating 
wavefunctions with identical periodicity as that of the SF phase, have the following 
form: 

V^,,(^,t)= /^cosV-* (5) 

y 92 



^,,(;^,t)= /-^sin^e— (6) 
\ 92 ^ 

with the corresponding energy i?/ = |^ \g\ — 3g2{2 + Vq — 2yu)]. In this case, the phase 
is arbitrary but the number of atoms on each lattice site is fixed. Thus, the atoms stop 
tunnelling between the lattice sites and the superfluidity is lost. The filling fraction, in 



^r_]_ _|_ cos(n7r)-l -| 

potential and shows modulations similar to the superfluid phase. 



this case, is found to he uj = ^[—1 ± — ]' which is independent of the lattice 



3. BEC with quintic nonlinearity 

For tight transverse confinement, strong interactions and low densities, the quintic 
nonlinear interaction term dominates the GP equation. In this case, we obtain only 
a superfluid phase with same periodicity of the lattice potential: ip{z,t) = i/oTe*'^^"*'^^*, 
where, cxi = (— I cos^| and cui = | — /i. Interestingly, this superfluid phase 
does not have any background term. Amplitude of the lattice potential and the three- 
body interaction strength need to be of opposite sign for the solution to exist. The 
supercurrent is found to be space-time independent: Ji = aiVi = Ci, where Ci = 
Here, vi is the superfluid velocity: vi = Xi = f^- The filling fraction is found to be a 
constant: 1^1 = ^ {-Vo/g2)'^^'^ ■ 



4. BEC with cubic nonlinearity 



The GP equation, in presence of pure two-body interactions {gi = g) is known to exhibit 
trigonometric solutions of the form |17j . 



i,{z,t) 



\ 



Vo 
9 



cos^ U e 



iX2{z)-ioJ2t 



(7) 



where the background a is a free parameter. The nontrivial phase, X2{z) 

0J2 -- 



C2 tan-i[^/l-^tan( 



/a(a-^) 



^ — fi + ga and the integration constant 
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{i-2fi-2a}){i-2n-2uj+2Va) ^ rpj^^ average atoms number density per lattice site is found to be 
related with the background a: psf = Ci — Here the filling fraction is independent 
of number of lattice site n. The non-trivial phase (c2 7^ 0) indicates the presence of 
superfluid component, with a constant flow density J2 = C2. The potential, density 
and the phase have the same periodicity in this case, which is very different from the 
previous cases, with two- and three-body interactions. The ground state energy of the 
superfluid state 

Tig (Wi SaVo 




- TT/x ( 2a - — ) + ttH ( a - ^ 1 , (8) 




9 

shows a branch cut aX, a — — Q and 02 = where air and a2T are the critical 
values of background. For ayr = 0, a;2 -|- // = | and at a2T = — , ^2 + A* = | + clq] 
at these values the super-current {J2 — C2) vanishes and superfluid phase transits to 
the insulating phase. The energy is continuous; however, the first derivative at the 
transition point (T), is discontinuous like the previous case, hence, indicating that 
the phase transition, is of first order. The wave function in the insulating phase at 

air is, i^hiz^t) = cos ze'''^^-^^^ with the energy, Ej = - 3Vog - 4). 

Corresponding number density is ^} = indicating that in this insulating phase, the 
interaction and potential strength have opposite signature. For the second insulating 

phase at a2T, ^1^12(^,1) = ^(^) sin ze-'^^-z^+^a)*, with, E] = ^{5Vog - 8/x + 4), and 

the average number density ^] = 

The interaction and potential strength, here, should have the same signature. For 
a positive Vq, these two phases respectively belong to attractive and repulsive regimes. 
It is worth emphasizing that for a fixed Vq and coupling both these insulating phases 
only exist for a given number density of atoms. In this well as in the previous 

cases, interestingly we find that the dispersion affects the density marginally, where it 
adds a constant term in the phase. 

Keeping in mind the possibilities of charge density type of ground states in one 
dimension, it is natural to investigate if other type of solutions, still commensurate with 
lattice periodicity, are allowed in this system. We have found such an exact solution 
exhibiting rational character: 

^,(,,,) = £±^£!^l!)±£^,-.., (9) 

1 -I- acos^ U 
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with a = (3 = 1 and 6 = 2. These solutions do not exist in superfluid phase. The 
parameters for this insulating phase are given by, 

3 rii 1 /12 91^ , PPi 

with, Go's = I — /i. Here, the competition of the lattice potential with the non-linearity 
is responsible for the non-perturbative rational solutions with dual frequency character, 
both commensurate with the lattice potential. It is worth observing that, in case of 
solitons, non-linearity and dispersion compensate each other, leading to stable localized 
solutions, as well as periodic non-sinusoidal cnoidal waves. Presently, the dispersion 
affects the character of the solution marginally. For non-singular solutions, Vq < — | in 
repulsive interaction regime {g > 0), and Vq > | in attractive regime {g < 0). It should 
be noted that, in contrast to the earlier found insulating phases, which are independent 
of the background, the present density wave solution has a background, which depends 
both on potential and the coupling. We further observe that, both the modulations are 
necessarily present in the presence of the lattice potential. 

The energy function is non-analytic as c? — >■ 0, ±1. For d = 0, the Fade type of 
solutions representing insulating phase reduces to the insulating phase exhibiting single 
frequency behaviour found earlier. In the limit d = ±1, the BEC is pushed gradually to 
the unstable position on top of the hills of the lattice potential indicating the unphysical 
nature of this singular point. 



5. Experimental Realization 

In the superfluid phase, the phase coherence between the different sites of the optical 
lattice in a BEC has been established experimentally through the formation of the 
interference pattern between condensates located at the nodes of the laser standing 
wave. In the insulating phase, the phase coherence is lost and hence, no interference 
pattern is formed. The phases identified here, can be detected in a method used by 
Cataliotti et. ai, in [15] for pure two-body interaction. The density wave insulating 
phase will not show phase coherence, however, unlike the regular insulating phase, the 
spatial distribution of atoms will show periodic behaviour commensurate with the lattice 
potential. The stripe phase will show phase coherence between the atoms, where the 
atoms will occupy only the alternating lattice sites. 

To check the parameter regime for the possible existence of these type of density 
wave and stripe phase, we consider the ^^Rb condensate with 10^ number of atoms. Mass 
of the atom is m = 1.44 x 10~^^ kg and the transverse trapping frequency is ujj_ = 2-7? x 140 
Hz. The s-wave scattering length a = 5.4 nm [26l [27] . The two-body interaction 
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strength is given by, gi = 2amu±/hk. The wave vector k = 2n/\ = 8.06 x 10^ m~^. It is 
worth mentioning that aU the above parameters can be changed independently. Larger 
values of interaction strengths can be achieved by higher density of the condensate, 
higher a and for small k. The scattering length can be tuned as well by Feshbach 
resonance. Small k may be achieved by adjusting the relative angle between the two 
interfering laser beams. The optical lattice depth is taken asVo — 0.002 (in the units of 
Er), implying that the potential is shallow, where the mean field GP equation is valid. 
As we increase Vq, at some point, phase coherence is lost and superfluidity breaks down 
and there will be no tunnelling of atoms between consecutive lattice sites. Hence, the 
phenomena discussed in this paper should be observable within the present experimental 
capability. 

6. Conclusion 

In summary, a superfluid density wave, having a periodicity twice that of the lattice 
potential, exists in the presence of both cubic and quintic nonlinear interactions, where 
the density modulation is over a constant background. A stripe phase is present in 
certain parameter regime, where superfluid matter is present in only some domains of 
the lattice. These phases arc absent, when, only cubic and quintic nonlinear interactions 
arc present. This density wave phase transits to an insulating phase through DSIT. The 
phase coherence is lost due to the modulational instability. The stability analysis of 
this solutions under VK criterion, supports the development of modulational instability 
in this system. The filling fraction is found to vary inversely with the lattice location 
from a given site, and differs for odd and even n. This different behaviour of the filling 
fraction in the neighbouring sites can be traced to the presence of the both cubic and 
quintic nonlinearities and is absent when pure cubic or quintic nonlinear interactions are 
operative. In pure two-body case, the periodicity of the density is same as the lattice 
potential. A density wave insulating phase is observed in this case, having two frequency 
modulations, commensurate with the lattice potential. Explicit calculations show that 
these type of density wave states and phase transitions can be probed in the present 
laboratory conditions. 
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